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ABSTRACT 



Context. The high sensitivity of the XMM-Newton instrumentation offers the opportunity to study faint and extended sources in the 
Milky Way and nearby galaxies such as the Large Magellanic Cloud (LMC) in detail. The ROSAT PSPC survey of the LMC has 
revealed more than 700 X-ray sources, among which there are 46 supernova remnants (SNRs) and candidates. 

Aims. We have observed the field around one of the most promising SNR candidates in the ROSAT PSPC catalogue, labelled 
[HP99] 456 with XMM-Newton, to determine its nature. 

Methods. We investigated the XMM-Newton data along with new radio-continuum, near infrared and optical data. In particular, spec- 
tral and morphological studies of the X-ray and radio data were performed. 

Results. The X-ray images obtained in different energy bands reveal two different structures. Below 1.0 keV the X-ray emission shows 
the shell-like morphology of an SNR with a diameter of ~ 73 pc, one of the largest known in the LMC. For its thermal spectrum we 
estimate an electron temperature of (0.49 ± 0.12) keV assuming non-equilibrium ionisation. The X-ray images above 1.0 keV reveal 
a less extended source within the SNR emission, located 1' west of the centre of the SNR and coincident with bright point sources 
detected in radio-continuum. This hard component has an extent of 0.9' (i.e. ~ 13 pc at a distance of ~ 50 kpc) and a non-thermal 
spectrum. The hard source coincides in position with the ROSAT source [HP99] 456 and shows an indication for substructure. 
Conclusions. We firmly identify a new SNR in the LMC with a shell-like morphology and a thermal spectrum. Assuming the SNR to 
be in the Sedov phase yields an age of ~ 23 kyr. We explore possible associations of the hard non-thermal emitting component with a 
pulsar wind nebula (PWN) or background active galactic nuclei (AGN). 

Key words, galaxies: Magellanic Clouds, ISM: supernova remnants 



1. Introduction 

The study of supernova remnants (SNRs) is crucial for a com- 
plete understanding of the chemical composition and evolution 
of the ISM in a galaxy because of their energy and matter inputs 
into the interstellar medium (ISM). The sample of SNRs stud- 
ied in our Galaxy is biased because of the high absorption in 
the Galactic disk. Therefore, observations of nearby galaxies al- 
low us to perform unbiased population studies and constrain the 
physical properties of the sources in detail. It is thus also possi- 
ble to better understand the evolution and structure of the ISM 
in the Milky Way. 

Located at a dista nc e of ~ 5 kpc to the Earth 
dFreedman et all l200lb iMacri et all [2006), the Large 
Magellanic Cloud (LMC) offers the ideal laboratory for 
studying a large sample of different types of objects (such 
as SNRs) in greater detail than in any other galaxy. Since 
its first detection in X-rays dMark et all I1969I) . the LMC has 
been extensively observed, but the major step forward came 
from more than 200 observations in a 10° x 10° field centred 
on the LMC, which have been performed with the ROSAT 
Position Sensitive Proportional Counter (PSPC) from 1990 



to 1994. For a description o f the ROSAT mission and PSPC 
detectors, see Triimper ( 1982J), iBriel & Pfeffermannl d!986l) and 
Pfefferma nn et all d 19871) . This survey revealed 758 sources 
(lHaberl & PietschL 1 1 999t hereafter labelled HP99), among 
which 46 sources were classified as firmly identified SNRs or 
candidates. 

Several SNRs in the LMC have been invest i gated using ob- 
serva tions with ROSAT (Filipovic et al., 1998; Williams et al 
1999 ), Chandra X -ray Observatory (iHughes et all 1 2006 
and XMM-Newt o n dWilliams et all [2004, 
ICrawford etail 120101) satellites, allow- 



Seward et al 



Klim ek et al 



* Based on observations with XMM-Newton, an ESA Science 
Mission with instruments and contributions directly funded by ESA 
Member states and the USA (NASA). 



ing a more d e tailed view of their morphologies and spectra. 
Badenes et al ] (120101) have studied the size distribution of the 
SNRs in the Magellanic Clouds (MCs), which has a maxi- 
mum at ~40 pc and may extend up to sizes of ~100 pc. With 
an extent of over 100 pc, SNR 0450-70.9 and SNR 0506- 
6542 (DE M L72) are among the largest SNRs detected in 
the LMC dWilliams et all 120041 ICaiko et all l2009HDesai et all 
2010; Klimek et all 1201 Oh . which may be highly evolved (age 
up to 100 kyr). The size distribution of the MC SNRs as well 
as those in our Galaxy or the neaby spiral galaxy M33 can- 
not be explained only by the Sedov expansion model for SNRs, 
but seems to be largely affected by the ambien t ISM densities 
dBandiera & Petruki l2010t iBadenes et all 12010, and references 
therein). 
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Furthermore, multi-frequency observations of several SNR 
candidates located in the LMC have enabled their f irm identifica- 
tion based on morpho l ogical and spectral criteria dBoiicic et al.L 
l2007t ICrawford et all [2010) and have revealed a strong cor- 
relation between the X-ray sources and the emission observed 
by the Magellan ic Cloud Emission Line Survey (MCELS; 
Smit h et all l2000). Indeed, an enhancement of the [S n] and Ho- 
coincident with the X-ray emission can be observed in most 
cases. In particular, the ratio fSnl/Hg is oft en higher than 0.4 
dLevenson et all IT995fc IWilliams et al., 2004). However, several 
SNRs in the Magellanic Clouds such as SNR LMC J0528-6714 
or/and SMC SNR JO 10505-7223 19 do not have any optical emis- 
sion. 

While the emission of the gas shocked by the shock waves 
of SNRs is mainly of thermal nature, there can also be a pul- 
sar or a pulsar wind nebula (PWN) in a SNR, which produce 
non-thermal emission. Pulsars are rapidly rotating neutron stars 
characterised by shor t periods (up to a few 10 s) and high sur- 
face magnetic fields dManchester et all 120051) . The dissipation 
of the rotational energy of pul sars via magnetised part icle winds 
can be at the origin of PWNe dGaensler & S land 120061) . A high 
percentage of PWNe known in our Galaxy have been detected 
in X-rays, and present a po wer-law spectra with a mean spec - 
tral index of T = -1.8 ± 0.6 dKargaltsev & PavlovL 120081 1201 Oh . 
Sensitive X-ray observations have enabled the detection of sev- 
eral PWNe and candidates within the Magellan ic Clouds, with 
similar properties as PWNe in the Milky Way dGaensler et all 
l200l IWilliams et all 120051 lOwen et al.Ll20lTh . 

Th ere are now o ver 5 well-establis h ed SNRs in the 
LMC dBadenes et all 12010b iKlimek et all l2010l and refer- 
ences therein) and som e additional ~ 20 SNR candidates 
(Boz zetto et allliii prep J) . This would comprise one of the most 
complete samples of SNRs in external galaxies. Therefore, it 
is of prime interest to study LMC SNRs and compare them 
with SNRs in othe r galax ies such as M 33 dLong et all l2010h . 
M83 dDopitaetalll2010h. the Small Magellanic Cloud (SMC ; 
IFilipovic et a ll 120051; iPavne et all [20071: iFilipovic et all 120081) 
laxydSti 



and our Galaxy ( Stu par et alll2008HGreenL l2009). 



The SNR candidates in the ROSAT PSPC catalogue 
have been classified based on the X-ray spectrum and 
spatial extent. Additional comparison to radio data taken 
with the Molo nglo Observatory S ynthesis Telescope (MOST; 
v=843 MHz; iTurtle & Miilsl 1 19841) and with optical data of 
MCELS has shown that there are ROSAT PSPC sources with ra- 
dio or optical counterparts indicative of an SNR but with a hard 
X-ray spectrum and therefore no typical SNR characteristics in 
X-rays. The source [HP99] 456 is one of the most promising 
new SNR candidates of this kind, with hard emission detected 
by ROSAT PSPC and a possible radio or optical counterpart. 
It has been recently re-investigated through new observations 
with XM M-Newton and the A ustralia Telescope Compact Array 
(ATCA; iHughes et alll2007h . 

In this paper, we report on the results of the analysis of new 
XMM-Newton and ATCA follow-up observations of the source 
[HP99] 456. Section [2] presents the observations and analysis 
techniques. Results of the multi-frequency analyses are pre- 
sented in Section [3] Section [4] presents the discussion on the 
SNR properties and on the different scenarios to explain the non- 
thermal emission. In particular, we explore a possible association 
of this second component with a PWN or an AGN. Conclusions 
are presented in Section[5] 



Fig. 1. Smoothed RGB three-colour image of combined 
exposure-corrected XMM-Newton EPIC pn, MOS1 and MOS2 
data (red: 0.2 - 1.0 keV; green : 1.0 - 2.0 keV; blue : 2.0 - 4.5 
keV, square root scale). The instrumental background was es- 
timated using the method described in Section 13.2.11 and sub- 
tracted from the images. The renormalisation factor was derived 
from the shaded detector corners. The field of the radio image 
presented in Figure[5]is overlaid for comparison (green square). 
The position of the ROSAT source [HP99] 456 is represented by 
a yellow circle. 



2. Observations and data reduction 

2.1. X-rays 

The XMM-Newton satellite is an X-ray observatory operated by 
the European Space Agency (ESA). The source [HP99] 456 has 
been proposed for observations with XMM-Newton (Obs. Id. 
0651880101, P.I. : M. Sasaki). This paper presents results of the 
X-ray analysis of a 20 ks observation obtained on this source. 

The source [HP99] 456 has been observed on 2010 June 
06 (from 04:18:56 to 09:52:50 UT) with the European Photon 
Imaging Camera (EPIC) in full-frame mode and thin filters. 
Using the EPIC MOS1, MOS2 and pn CCDs, it offers the op- 
portunity to perform sensitive X-ray observations of a field of 
the sky of diameter of 30'. More detail ed technical descrip tions 
of the EPIC camera s are presented by iTurner et al.l (12001 1) and 
IStriideret al.ld2001l) . 

The EPIC data were analysed with SAS vl0.0.(fl The expo- 
sure time, after removing periods of high background, is ~ 1 1 ks. 
Pixels flagged as bad were not taken into account and screening 
on the patterns (from to 12 for MOS; from to 4 for pn), corre- 
sponding to the canonical set of valid X-ray events (calibrated on 
the ground) was applied for the image and spectral analyses. We 
performed source detection by using the SAS tasks eboxdetect 
and emldetect. 



Science Analysis Software: http://xmm.esac.esa.int/sas/. 
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2.2. Radio-continuum 

The field of [HP99] 456 was observed with the Australia 
Telescope Compact Array (ATCA) on 2010 November 29, with 
an array configuration 6C, at wavelengths of 3 and 6 cm (9000 
and 5500 MHz), and a bandwidth of 2 GHz (ATCA project 
C2367). The observations were carried out in snap-shot mode, 
totalling about 1 hr of integration over a 12 hr period. PKS 
B1934-638 was used for flux and bandpass calibration and PKS 
0530-727 was used for phase cali bration. Standar d calib ration, 
editing and imaging techniques (Sault & Kille enj | 201 lb were 
used. Large bandwidth multifrequency clean (Sault & Wieringa, 
1994) was used to deconvolve the image. We point out that inter- 
ferometers such as the ATCA suffer from missing flux owing to 
the lack of short spacings, which significantly affects the overall 
detection of extended emission like that from SNRs. 

We also used various other radio observ ations (See Table [TJ 
including 84 3 MHz by iMills et all lfl98l . 1377 MHz from 
iHughes et all d2007l) and 4800 MH z & 8640 MHz taken from 
the mosaic project of lDickell (|2005). 

2.3. Optical 

The Magellanic Cloud Emission Line Survey (MCELSfl was 
carried out from the 0.6 m University of Michigan/CTIO 
Curtis Schmidt telescope, equipped with a SITE 2048 x 2048 
CCD, which gave a field of 1.35° at a scale of 2.4"pixeL' 
dSmith et all l2006t IWinkler et ail lin preph . Both the LMC and 
SMC were mapped in narrow bands corresponding to Ha, [O m] 
(,1=5007 A), and [Sn] (^1=6716,6731 A), plus matched red and 
green continuum bands. All data were continuum subtracted, 
flux-calibrated and assembled into mosaic images (a small sec- 
tion of which is shown in Figure|6]l. 

3. Data analysis 

3.1. X-ray morphology 

The morphological analysis was performed in three energy 
ranges : 0.2 -1.0 keV (soft band), 1.0-2.0 ke V (medium band), 
2.0 - 4.5 keV (hard band). For each energy interval, exposure- 
corrected images were obtained using SAS tools using the fol- 
lowing method. 

First, X-ray data were cleaned using the selection procedure 
described in the previous section. We computed event images 
and exposure maps using the filtered dataset. Then we applied 
a mask to remove bad pixels from the three instruments. We di- 
vided each resulting event image by the corresponding exposure 
map and smoothed the resulting images with a Gaussian filter 
using asmooth. Finally, we added the event images from the dif- 
ferent instruments for each energy band. 

Figure Q] shows the RGB exposure-corrected image around 
the unidentified source [HP99] 456 obtained with the EPIC data 
in the three bands defined above (red: 0.2 - 1.0 keV; green : 
1.0 - 2.0 keV; blue : 2.0 - 4.5 keV) after subtraction of the in- 
strumental background (see caption). The emission seen in the 
hardest band is similar to the medium band. This morphological 
analysis reveals the existence of two distinct emitting regions: 

1. A soft emitting region: below 1.0 keV, the X-ray emis- 
sion is dominated by a bright thermal (for more details, see 
Section [3~2b component. It presents a shell-like morphology 
with an extent of ~ 5.0'x 4.6'. 

2 MCELS : http://www.ctio.noao.edu/mcels/ 



2. A hard emitting region: above 1 .0 ke V, the X-ray analysis re- 
veals a less extended harder component. Its centre is located 
~ 1' west from the centre of the soft shell-like emission and 
presents an elongated morphology of length ~ 0.9'. The po- 
sition is consistent with the ROSAT position of [HP99] 456. 

The presence of the hard emitting region coincident with the 
ROSAT source suggests that [HP99] 456 is unrelated to the shell- 
like emission region. Therefore, we designated the soft, shell- 
like source XMMU J0541. 8-6659 according to its approximate 
central coordinates. Additional details on the positions and ex- 
tents of both emitting regions are summarised in Table [2] The 
count rates of the EPIC-pn observations are listed in the last col- 
umn. 




Fig. 2. Smoothed image of exposure-corrected XMM-Newton 
EPIC data in the 0.2 - 1.0 keV energy range. The dashed cir- 
cle shows the region considered for the estimation of the local 
background. The thick circle represents the region used to de- 
rive the spectrum of the soft emission from the shell-like struc- 
ture. Cancelled circles show the position of the bright sources 
excluded from the spectral analysis. 



3.2. X-ray spectral analysis 

The following sections present the results of the subsequent 
spectral analyses of the two spectral components mentioned 
above. We first describe the method used to estimate the contri- 
bution from the intrinsic detector and X-ray background for each 
emitting region. The spectral analys i s was performed using the 
XSPE C v 12.6.0 package dArnaudl ll996UDorman & Arnaudl 
I2001h . 

Point sources that were detected using the method described 
in Section lXTI were excluded from the spectral analysis. We used 
data between 0.3 keV and 5.0 keV because no significant emis- 
sion is detected at higher energies. 
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Table 1. Integrated radio flux densities of [HP99] 456. 



V 


A 


R.M.S 


Beam Size 


STotal 


Reference 


(MHz) 


(cm) 


(mJy) 


(") 


(mJy) 






843 


36 


1.5 


43.0x43.0 


0.045 


Mills et al. (1984) 


1377 


20 


1.5 


45.0x45.0 


0.029 


Hughes et al. (2007) 


4800 


6 


0.6 


30.0x30.0 


0.008 


Dickel (2005) 


8640 


3 


0.6 


12.0x12.0 


0.004 


Dickef (2005) 


5500 


6 


0.1 


2.8x2.2 




This Work 


9000 


3 


0.1 


2.8x2.2 




This Work 



Table 2. Morphological details on the two emitting regions. The last column lists the count rates estimated from the analysis of 
EPIC-pn data between 0.3 and 5.0 keV. 



Emitting 


Ra 


Dec 


Size 


Position angle 


pn count rate 


region 


(hh:mm:ss.d) 


(V:") 


(') 


o 


(cts s~') 


Soft emission 


05:41:51.5 


-66:59:02.8 


5.0x4.6 


-45 


(18 ± 2) x 10~ 2 


Hard emission 


05:41:39.4 


-66:58:45.8 


0.9 x 0.45 


-15 


(2.8 ± 0.3) x 10" 2 



3.2.1. Estimation of the instrumental background 

We estimated the contribution from the intrinsic instrumen- 
tal background using the filter wheel closed (FWC) data 
(Freyberg eTaILl2001albh . The intrinsic background is composed 
of internal electronic noise as well as continuous and fluorescent 
X-ray emission induced by high-energy particles. This back- 
ground is measured by operating the EPIC instruments in the 
"closed" filter wheel position, where no photons from astrophys- 
ical sources can be observed. The FWC spectrum needs to be 
renormalised using the continuum at higher energies. 



3.2.2. Estimation of the X-ray background 

The contribution of the X-ray local background emission was 
derived from the same observation as the source itself. We de- 
fined a region located close to our source (see Figure 0, and 
extracted the spectrum from the observation and the FWC data. 
We estimated the renormalisation factor of the FWC spectrum 
with respect to the observed spectrum using the total count rates 
measured above 5 keV in the corners outside the field of view of 
the EPIC camera. The value of the renormalisation factor is 1.3. 

The spectrum of the diffuse emission after subtracting the 
instrumental background can be modelled as a sum of the fol- 
lowing three components: 

- a soft component, modelling the emission from the local 
bubble with a collisional plasma, non-equilibrium model 
(NEI: iBorkowski et all 11994 1200 It lHamilton et allll983b 
iLiedahl et al.L 1995 J3 assuming a low temperature (0. 1 keV); 

- a hard component modelled with an absorbed NEI model 
with a temperature of 0.3 keV, which stands for the emission 
from the Galactic halo; 

- an absorbed power-law of spectral index F = 1.46, to de- 
scribe the extragalactic unresolved emission. 

The parameters are fixe d to values similar to those used by 
iKuntz & Snowdenl(l2010l) . 

The two last components are convolved with a Tiibingen- 
Boulder ISM absorption model (TBABS). The Galactic fore- 
ground hydrogen column density of 0.6 x 10 21 cirT 2 was derived 



from Star k et all (1 19921) and was used for the absorption in the 
Galaxy. 

3.2.3. Soft emission region 

The X-ray emission below 1.0 keV is dominated by a struc- 
ture presenting a shell-like morphology with a maximal di- 
ameter of ~ 5.0' and centred at RA(J2000) = 05 ;, 41 m 51.5*, 
DEC(J2000) = -66°59'03.8". Photons located within the region 
corresponding to the hard X-ray emitting region, as defined in 
Table [2] were excluded from the source spectrum to avoid any 
contamination, as were detected point sources visible in Fig.Q] 

The spectrum of the shell-like structure was obtained after 
the subtraction of the instrumental background. The X-ray back- 
ground was estimated using a local background spectrum ex- 
tracted from the same data. The source spectrum and the local 
background were modelled simultaneously in XSPEC by using 
a background model consisting of the components as explained 
in the previous section. All parameters of the three background 
components were fixed except the normalisations. The normali- 
sations of the background components and the parameters of an 
additional model component for the shell emission were free fit 
parameters. 



Table 3. XMM-Newton spectral results of the soft thermal emit- 
ting region in the 0.3 - 5.0 keV energy range. Model is VNEI 
with abundances fixed to 0.5 of the solar abundances. 



Parameter 


Value 


N H (xlO 21 cnT 2 ) in the Galaxy 


0.6 (fixed) 


N„ (xlO 21 cm" 2 ) in the LMC 


(fixed) 


kT (keV) 


0.49 ± 0.12 


Ionisation timescale* t (s cirT 3 ) 


(1.5+0.7) x 10 10 


Absorbed flux (0.3 - 5.0 keV, erg/cm 2 /s) 


2.5 x 10~ 13 


Absorbed luminosity (0.3 - 5.0 keV, erg/s) 


5.1 x 10 34 


Unabsorbed luminosity (0.3 - 5.0 keV, erg/s) 


8.9 x 10 34 


Reduced^- 2 


1.27 


Degrees of freedom 


86 



3 For more details on XSPEC models, please 
http://heasarc.nasa.gov/xanadu/xspec/manual/manual.html 



*: The ionisation timescale r = n e t, where n f is the electron number 
density and t is the age of the gas. 
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The spectral analysis of the soft circular emitting re- 
gion reveals a thermal spectrum. It can be modelled 
with a single-temperature non- equilibrium ion i sation colli- 
sional plasma model (VNEI; j orkowski et all 11994 l200lt 
iHamilton etall 1 1 9 8 3b iLiedahl et allll995l) with "temperature of 
(0.49 + 0.12) keV with an absorption of N H = 0.6 x 10 21 cnr 2 
for the Galactic foreground (with solar abundances) and a value 
of Nh for the absorption in the LMC. The abundances were fixed 
to 0.5 of the solar system abundances for the emission and ab- 
sorption taking place in t he LMC. This value is t he standard 
mean value for the LMC dRussell & Dopital 119921) . The value 
of Nh for the LMC turned out to tend towards with an 90 % 
C.L. upper limit of 0.4 x 10 21 crrT 2 . This parameter was there- 
fore fixed to 0. 

The resulting spectrum is presented in Fig. [3] The corre- 
sponding spectral parameters are given in Table [3] 



I = /d.D.f. = 1.27(B6d.0.1.) 




0.5 1 2 

Energy (kBV) 

Fig. 3. XMM-Newton EPIC-pn spectrum and model of the soft 
emitting region with a VNEI model fit (black points and solid 
line). The estimation of the local background spectrum and 
model (red points and solid line) is described in Section I3.2.1I 
The background spectrum was not subtracted from the source 
spectrum but was modelled simultaneously and is included in 
the spectral model of the source spectrum. The different spectral 
components are shown separately (red and black dotted lines) 
for the background and source spectra. The bottom panel shows 
the residuals from the best-fit models. 



3.2.4. Hard emission region 

The X-ray emission above 1.0 keV is dominated by a narrow 
structure of length ~ 0.9' located ~ 1' to the west of the centroid 
of the shell-like thermal emission, as seen from Figure[T] Above 
1.0 keV, the hard source has a surface brightness of ~ 5.8 x 10~ 3 
counts/s/arcmin 2 , i.e. ~ 6 times the background level, which has 
a surface brightness of ~ 9.4 X 10~ 4 counts/s/arcmin 2 . There is 
indication for substructure in the X-ray image, which could be 
caused by multiple sources. The ROSAT source [HP99] 456 co- 
incides with the northern part of the narrow structure. 

To analyse the spectrum of this second component again the 
instrumental background spectrum was subtracted and the X-ray 
spectrum was modelled simultaneously, this time using the sur- 



I ■ I d.o.1. = 0.22 (6 d.o.f.) 
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Fig. 4. XMM-Newton EPIC-pn spectrum of the hard emitting 
region with a power-law model fit. The solid line is for the hard 
emission plus background spectra. Here, the background is the 
sum of the shell emission and the local emission. The back- 
ground spectrum and shell spectrum were evaluated as described 
in Section I3.2.4I The background spectrum was not subtracted 
from the source spectrum but was modelled simultaneously and 
is included in the spectral model of the source spectrum. The 
contribution from the hard emitting region is represented in red. 
The bottom panel shows the residuals from the best-fit model. 



rounding shell-like emission as the local background. The nor- 
malisations of the component corresponding to the X-ray back- 
ground components were renormalised to the area and were fit- 
ted, along with the spectral parameters of the hard emitting re- 
gion. 

The spectrum of the hard emitting region is non-thermal and 
can be modelled in the 0.3 - 5.0 keV energy range with an ab- 
sorbed power-law with a spectral index of F = 1.8+0.3, with an 
absorption of Nh = 0.6 x 10 21 cm~ 2 for the Galactic foreground 
(fixed, with solar abundances) and a value of Nh for the absorp- 
tion in the LMC of N H = (0.58 ± 0.53) x 10 21 cm 2 . This yields 
an absorbed flux of ~ 9.5 x 10~ 14 erg/cm 2 /s (reduced^ 2 = 0.22 
for 6 degrees of freedom). The resulting spectrum is shown in 
Fig- El In particular, the contribution from the hard X-ray emit- 
ting region is represented in red. 

An alternative technique was used to derive the spectrum of 
the hard emitting region. Instead of using the FWC data, we 
extracted a background spectrum from a nearby region located 
within the shell and subtracted it from the source spectrum. This 
method yields consistent results for the hard emission. 

3.3. Radio 

New high-resolution (~2") radio-continuum observations re- 
solved low-resolution emission into three bright sources within 
the hard X-ray emitting region (Fig. [5]). For comparison, the size 
of the radio image shown in Fig.|5]is represented in Fig.Q] In par- 
ticular, three bright sources can be seen coinciding with the cen- 
tral position of this hard non-thermal X-ray emitting region. The 
brightest source (the most eastern one) has a spectral index of 
or = -0.5 + 0. 1 with integrated flux densities of (0.85+0.05) mJy 
at 5.5 GHz and (0.67+0.05) mJy at 9.0 GHz. 
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T| ■0.350 




5 h 41 m 45 s 40 s 35 s 30 s 

Right Ascension (J2000) 

Fig. 5. High-resolution ATCA radio-continuum image at 6 cm 
(5.5 GHz) overlaid with 20 cm contours (blue) from the low- 
resolution mosaic image of the immediate surroundings of 
[HP99] 456. Contours are from 3 to 18 mJy/beam in steps of 
1 mJy/beam and 6 cm beam size (2.8" x2.2") is shown in the 
bottom left corner. 



Our estimate of the SNR overall radio spectral index (exclud- 
ing central point source; using ux densities from mosaic surveys 
only and listed in TableQ]i.e. observations from 843, 1377, 4800, 
8640 MHz) is a — -1.0 ± 0.2, which is more typical for younger 
SNRs. 

3.4. Optical 

The emission from optical lines (He and [S n]) derived from the 
MCELS were considered. The Ha emission in the surroundings 
of [HP99] 456 are presented in Fig. [6] This image shows a pos- 
sible correlation between the western part of the shell-like X-ray 
thermal emission observed by XMM-Newton, the low -resolution 
radio-continuum and the Ha emission line, as seen in Figures 
[5] and [6] However, we acknowledge that this correlation could 
be a chance coincidence because the Ha "cloud" could be an 
unrelated Hn region. 

Low-resolution radio-continuum images at 36, 20, 6 and 
3 cm show good alignment with the optical (MCELS) feature 
and also coincide with the maximum emission from the X-rays. 

Supernova remnants can show significant Ha, [Sn] and 
[Om] line emission if they are in the radiative phase. In this 
case, the flux ratio [S n]/Ha is a crucial parameter to dis ti nguis h 
between SNRs and e.g., Hn regions dLevenson et all Q995). 
However, the [S n]/Ha ratio around [HP99] 456 is ~ 0.3, which 
is not necessarily indicative of SNR emission. 

Ha emission is spatially coincident with the hard emitting 
region, but additional studies are required to determine if they 
are associated or not. 

4. Discussion 

Thanks to the high sensitivity of XMM-Newton, a recent observa- 
tion of the field of [HP99] 456 has revealed the existence of two 
emitting components. The softer component presents a circular 




Fig. 6. Ha emission (arbitrary units) derived from the MCELS 
in the surroundings of [HP99] 456. A contour of the soft ther- 
mal emission of XMMU J0541. 8-6659 is overlaid in green for 
comparison. 




Fig. 7. Soft X-ray emission of XMMU J0541. 8-6659. The con- 
tours of the near-infrared emission (Spitzer, 5.8 //m) are overlaid 
for comparison. 



morphology and is firmly identified as a new SNR, designated 
XMMU J0541. 8-6659. The nature of the emitting region with 
harder spectrum is still unclear, but we examined the possible 
identification as a PWN or a background AGN. 
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4.1. Properties of the newly identified SNR 
XMMU J0541. 8-6659 

X-ray data reveal a large structure of soft emission close to the 
ROSAT PSPC source [HP99] 456 in the LMC. Analyses per- 
formed in different energy ranges indicate that this source dom- 
inates the X-ray emission below 1.0 keV. 

The spectrum of the SNR XMMU J0541 .8-6659 can be mod- 
elled with an absorbed component representing a plasma in non- 
equilibrium ionisation characterised by a temperature of ~ 0.49 
keV, which is the average of the typical temperature of previ- 
ously detected SNRs in the LMC. We note that the spectral pa- 
rameters quoted in Table [3] are consistent with the parameters of 
the unambiguously identified SNRs in the LMC dKlimek et all 
lIoTotlLevenson et all [T99l I Williams etafl [1999112004 . It fur- 
ther supports the identification of this thermal emission with a 
SNR. 

The soft emitting region is firmly identified as a new SNR, 
in view of its morphology and X-ray spectrum. It has a max- 
imum diameter of ~ 5.0', which corresponds to a maximum 
extent of ~ 73 pc at a distance of ~ 50 kpc. Therefore, 
XMM UJ0541.8-6659 is one of the largest SNR observed in the 
LMC (IWilhams et all [19991: iKlimek et alllloToh . 

However, the large extent of the source does not necessarily 
imply an old age for the SNR. Indeed, assuming the SNR to be in 
the Sedov phase, the dynamic age of the source can be estimated 
using the shock temperature as follows : 

k B T s = 1.8 x 10 5 (y) keV, (1) 

where k B is the Boltzmann constant, T s is the temperature of 
the shock, which is comparable here to the plasma temperature 
estimated in Section 13.2.31 R s is the radius of the shock (in pc) 
and t is the age of the SNR (in yr). 

From the spectral fit yielding a plasma temperature of k B T s ~ 
0.49 keV, we derived an age of ~ 23 kyr. This value is well below 
the a ge of previous large l y-extended SNRs dete cted in the LMC 
(e.g. lKlimek et alll2010tlWilliams et alll2004b . 

This result means that the large extent of the source is not due 
to its high age, but rather to an expansion in a probably rarefied 
ambient medium. This is not surprising knowing the position of 
XMMU J0541. 8-6659 at the north-east part of the LMC, where 
the gas density is quite low. 

As a legacy project of the Spitzer Space Telescope, a survey 
of the Magellanic Clouds has been carried out cal led "Surveying 
the A gents of Galaxy Evolution" (SAGEfl (Meixner et al., 
2006) to obtain images and spectra of the dust emission. A uni- 
form survey of the LMC was performed in a 7° x 7° by the 
Infrared Array Camera (IRAC, dedicated to near infrared obser- 
vations at wavelength 3.6, 4.5, 5.8, and 8 pm) and Multiband 
Imaging Photometer for Spitzer (MIPS, mid- and far-infrared 
observations at wavelengths 24, 70, and 160 pm). The near- 
infrared emission (at 5.8 pm) in the region surrounding the 
sourc e [HP991 456 is presen ted in Fig. [7] Spitzer-S AGE observa- 
tions ( Me ixner et all 120061) reveal dust emission spatially coin- 
cident with the eastern half of the shell-like soft emitting region. 
However, there is no clear indication that this dust emission is 
associated to the source XMMU J0541. 8-6659. 

4.2. Hard emission region: a PWN or background AGN? 

Energy-dependent morphological studies of the X-ray emission 
in the surroundings of [HP99] 456 have revealed an extended 

4 Spitzer-SAGE: http://sage.stsci.edu/index.php 



elongated structure located within the SNR, which dominates 
the X-ray emission above 1.0 keV. Its centre is located ~ 1' from 
the centre of the SNR. This hard emission region has an appar- 
ent extent of ~ 0.9', corresponding to ~ 13 pc at a distance of 
~ 50 kpc and a non-thermal spectrum. However, the angular res- 
olution of XMM-Newton does not allow us to clarify if this hard 
source is indeed extended (with substructure) or if it is a con- 
glomerate of point sources. A careful look at the ROSAT images 
reveals that the position of [HP99] 456 is consistent with the 
hard source detected with XMM-Newton and in particular with 
the more northern knot in the hard emission (see Fig. U}. Hence, 
a possible association between the soft and the hard components 
and the exact nature of the hard X-ray source remains unclear. 
The radio observations reveal several point sources at the po- 
sition of the northern X-ray knot. This may point towards an 
identification as a pulsar with PWN. However, equivalently, one 
may argue for an existence of several background sources such 
as AGN. We investigate the different scenarios below. 

Firstly, we examine the possibility that the hard X-ray emit- 
ting region is associated to background AGN. With a flux of 
5.1 x 10~ 14 erg crrT 2 s _1 (0.5 - 2 keV), we expect to see ~ 12 
AGN per square degree in the sky dRosati et all 12002). We have 
an 8 % probability for an AGN inside a 5.0'x4.5' SNR shell. This 
low probability and the non-detection of any similar morpholog- 
ical (jet-like) structure either in radio, IR or in optical renders 
an association of the hard X-ray component with an AGN quite 
unlikely. 

It is worth noticing that the different values of absorption ob- 
tained during the spectral fit of the hard and soft emitting regions 
may indicate that they are located at different distances from us. 
While the absorption of the shell is consistent with the SNR be- 
ing located in the LMC, the slightly higher absorption of the 
hard emission suggests that the source might be located behind 
the SNR. 

Secondly, we consider the possible association of the hard 
X-ray component with a PWN. Three radio sources are embed- 
ded in the hard and narrow X-ray component. The brightest radio 
source presents a spectrum that can be modelled with a power- 
law of index a » -1.0 which is consistent with spectra of pul- 
sars. If this radio source is indeed a pulsar, then it might power 
a PWN that would be seen in the X-ray domain. Timing obser- 
vations are required to confirm or invalidate this identification. 
Furthermore, it is worth noticing that with a spectral index of 
r = 1.8 + 0.3, the spectrum of the hard X-ray emitting source is 
similar to spectra of PWNe in the Milky Way (spectral i ndices 
ranging from -2.5 to -1.2; iKargaltsev & Pavlovl I2008L 1201 Ob 
but also typical for AGN. 

The multi-frequency observations along with the spectral 
and morphological results of our X-ray and radio data analyses 
do not uniquely identify the nature of the hard X-ray emission 
embedded in the SNR XMMU J0541. 8-6659. Deeper X-ray ob- 
servations with higher angular resolution are required to differ- 
entiate between the two scenarios. 



5. Conclusions 

A recent X-ray observation with XMM-Newton of [HP99] 456 
previously discovered by ROSAT led to the identification of a 
new SNR (XMMU J0541. 8-6659) in the LMC, which presents 
a shell-like morphology and a soft thermal spectrum. No clear 
correlation can be found between the X-ray and optical emission. 
The physical properties (temperature, size, etc.) of this new SNR 
are consistent with SNRs previously identified in the LMC by 
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X-ray observations dLevenson et all Il995t iKlimek et all l2010t 
I Williams et aij.ll999[l2004lT ~ 

Additional analyses of the XMM-Newton data have revealed 
a harder and narrower emitting region within the shell, which 
may be a PWN or background AGN. This emission is likely the 
counterpart of the hard source discovered by ROSAT. Follow-up 
observations with ATCA revealed several radio point sources 
coincident with the hard emitting region, one of which could be 
a pulsar. Deeper observations of the hard X-ray emitting region 
with the Chandra X-ray Observatory will help to distinguish 
between the different scenarios and to unveil its nature. 
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